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Abstract 

In this paper we use abstract vector spaces and their duals without any canonical basis. Some 
of our results can be extended to infinite dimensional vector spaces too, but here we consider only 
finite dimensional spaces. We focus on a general perturbation problem. Assume that B : V V 
is a linear operator, which is perturbated to B' = B + Q. We examine the question how the 
determinant and the inverse change, because of this perturbation. In our approach the operator 
Q is given as a sum of dyadic products Q = X^iLi ® Pi, where Vi € V and pi £ V*. In this 
paper we derive an m-th order (m £ N) approximation formula for det B' and {B')^^ , which gives 
the exact result if m > fc. 

Introduction 

In this paper we use abstract vector spaces and their duals without any canonical basis. Our 
notations follow the classical linear algebraic notations, for details see for example [T]. We assume 
that 1^ is a real or complex vector space and the maps between vector spaces are linear. Some of 
our results can be extended to infinite dimensional vector spaces too, but here we consider only finite 
dimensional spaces. 

We focus on a general perturbation problem. Assume that B : V ^ V is a linear operator, which 
is perturbated to B' = B + Q . We examine the question how the determinant and the inverse change, 
because of this perturbation. One natural approximation is given by the Taylor expansion, but that 
process requires norm on the vector space, and the Taylor polynomials do not give the exact result 
when the series is cut within a finite number of terms. In our approach the operator Q is given as a 
sum of dyadic products 

k 
i=l 

where Vi G V and pi e V* for every i — 1, . . . ,k; moreover we can assume that k < dim V^. In this 
paper we derive an m-th order (m £ N) approximation formula for dcti?' and {B')^^, which gives 
the exact result if m > fc. 
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1 On the inverse of perturbated operators 

If A : V X V* is a. bilinear map, then there exists a unique k{A) : V ^ V map such that 

p{k{A)x) ^ A{x,p) y{x,p)eVxV*. 
We define the map k : hm(V x y*,]R) Lin(y, V), which is an isomorphism. 

Definition 1. Let V be a vector space, such that dimF = n > 2. For given vectors zi,. . . ,Zk G V 
and covectors p\,...,pk G V* , where 1 < < n — 1 define the map 

^ll-Sl-V XV* {v,q)^{qAp,A---Apk){v,z,,...,Zk). (1) 

We have k G Lm{V, V), and we introduce the symbol 

k 

B^izi,pi) = {zuPi) H(02,P2) □ . . . Hzk,Pk) = K • (2) 

It is obvious from the definition that if or (pi)i=i,...,fc are not linearly independent, then 

fc 

B{zi,pi) = 

and if tt, tt' are permutations of the set {!,...,/;} then we have 

fe , k 



1=1 1=1 



Moreover if v £V and q gV* then the equalities 

q 



.H(2;i,Pi)^ (v) = (qApiA---ApkY'V^Zi^---'Zk) 



/ fc \ 

f H j (v) = v{pi A ■■■ Apk){zi,...,Zk) - ^ZiipiA--- Apk){zi, . . . , Zi-i,v, Zi+i, . . . , Zk) 

^ i=l 

hold. 

Now we show how one can compute the determinant and the inverse of the perturbated identity 
operator. 

Theorem 1. Let V be a finite dimensional vector space, dimF = n and (Mi)i=i,...,fe a family of 
vectors and (pi)i=i,...,fe a family of covectors. Define the linear map 



fc 



^ = idy (3) 



We have for the determinant of A 



miii(n,fc) 

detA=l+ (a P,, ) (4) 



i=l l<]i<j2<---<]i<k 

and if det A ^ 0, then we have for the inverse of A 



min(r?,— l,fc) 

A-UetA = idv+ Z (5) 



i=l l<h<h<--<ji<k 
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which can be written in the following form: ifx^V and g G then 



min(n— l.fc) 

q{A-^x)detA = qix)+ ^ ^ 



i=l l<ji<32<---<ji<k 



1 = 1 



(6) 



Proof. To prove Equation (jj]), we assume that the vectors are hnearly independent and we 

complete it with elements {ui)i=k+i,...,n to get a basis in V, and assume that lj G A"{V) is a nonzero 
n-form. We compute the determinant from the following equation 

Lu{Aui, . . . , Aun) — lli{ui, . . . , u„) dct A. 

A simple expansion of the expression 

k k 
W I Ul +y^MjPj(ui), . . . ,U„ + y^MjPj(M„) 

gives Equation (jH). 

Let us denote the right hand side of Equation ([SJ by ip, define m niin(n — 1, fc), and assume that 
rn > 1, since the m = I case is trivial. If a; = Ay then we have the following equation. 



V ^q{y) + ^q{ua)pa{y) + ^{q A Pj){Ay,Uj) 

a=l j=l 

m—1 



E E 

i=2 l<ji< - <3,<k 



qA [ A^p,, 



{Ay,Uj,,...,UjJ 



q A [ A^pi 



{Ay,ui,...,Uk) 



We can expand the third term 

'^{qApj){Ay,Uj) = q{y)J2Pjiuj) " E '^("j)^^^ + EE^^<=(2/)(9 Apj)(uc,^ 



the summand in the fourth term 



9 A I A^p,, 



(Ay, , . . . , J = q{y) 



qA [ A^pj, 



{Uj^,...,Uj,) + 



c=l 



9 A ( A^Pj, 



{Uj, , . . . , J 



and the fifth term 



qA [ A^pi 



E P'^^y') 

?{/iT-.-,3a 



{Ay, ui,...,Uk) =q{y) { A^pi ] (ui, . . . , Uk) 



c=l 



. . . 
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Combining these terms we get 



i=2 
k 



(lie, '^ji 1 • • • •) "^ji ) 



+ E E E pM 

1=1 i<ji<---<j.<fc <^=i . 

c?tji . --.Ji} 

m i 

+ E E E(-i)^p>(y) 

4=2 l<ji< - <Ji<fc c=l 



9 A I A^Pj, 



[Uj-^ , . . . , Uj; 



This can be rewritten as 
ip = q{y) detA 



m—1 ( k r / ■ \ 

E E E p^(y^ '^^(Ap^') 



+ E E(-i)^?'.c(2/) 

l<Jl<---<j", + l<fc c=l 



/ i+1 



(Uc, Mjj , • • • , Itji ) 



(Ujj , . . . , Uj^^j ) 



(7) 



For a given l<j<m — 1, we assume that 1 < ji < • • • < < fc are fixed. Now we check how the 
term 

9a(^A^P,,^ (8) 
occurs in the previous formula. From the first summation in Equation ([7]), we have 

A; 

-. I i \ 

(lie; '^jl 7 ■ ■ • ' ^ji ) 



9 A ( A^Pj, 



(9) 



E 

c=l 

Ji} 

If in the second summation we have the indices 1 < < ■ ■ • < < k, we get the term ^ if 

(ill • • ■ J Jb-l5 Jb+li • ■ ■ J Ji+l) — (ill ■ ■ • 5 ji) 

holds for a b index. This b is the extra j' index in the second summation. To get the term ([S]) we 
have c = &, so the second summation is 



E (-^)''PrM 

t>({jl Ji} 

This can be written as 

k 



gA( Ap,; 



(Wjj, . . . . --Uj' ). 



E p^o^y^ 



9 A ( A^Pj, 



(10) 



l>^{ii,...,ii} 



Adding Equations ()9|10p . we get zero, which means that the summands for all 1 < i < jti — 1 in 
Equation ([7]) are zero. This proves that ip = q{y)detA, which is the left hand side of the Equation 
dH). □ 
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We note, that the previous theorem can be proved by induction on fc, but the detailed proof has 
approximately the same length. Now we can state our main result as a simple consequence of the 
previous theorem. 

Corollary 1. Assume that B : V ^ V is an invertible map and consider the perturbated operator 

k 

S' = B + ^v, (11) 
where Vi G V and pi G V* for all i = 1, . . . ,k. Let us define Ui = B~^Vi (for all i = 1, . . . ,k) and 

k 

A = idv+^u^>S'P^. (12) 

1=1 

// det A =^ 0, then we have for the inverse of the perturbated operator 

/ min(n-l,fe) . \ 

<«''"^d^ni«" + 5j:4 S S U-v (13) 

\ i=l l<il<j2<-<J,<fe / 

where n — diniT^. 

Proof. Since B' — BA, we use the formula {B')^^ — A^^B^^, where A~^ is given by Equation 

®. □ 

2 Connection with the Taylor expansion 

In applications, we assume that the perturbation 

k 

B' = B + Y.v^®p, = B + Q (14) 

1=1 

is small in some sense with respect to B. (In our framework there is no norm, so the word small has 
just intuitive meaning here.) If we take into account only m (to G N) (or a less number of) products 
of {ui,pi), we get the m-th order approximation of (B')^^, that is 

1 1 / i \ 

\«=1 l<3i<j2<-<Ji<k j 

where 

det„A=l + ^ai, a; = ^ ( A pj, j (uj^, . . . , UjJ and Ui = B~^Vi. 

We have for the zeroth, the first and second order approximation of {B')^^ 

{BTi' = - -^B-^QB-^ 
1 + ai 

= - ^±^B-'QB-' + \ B-'QB-'QB-' 
1 + ai + q;2 1 + ai + a2 
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and in general 



{BX: = + ^ X Y^i^iy 1 + ^ a, {B-'QrB-\ (15) 

It is obvious from the construction that [B')^ = {B')~^ \i m>k. 
Since the z-th derivative of the inversion function i{B) = B~^ is 

d't(B)(QH) = {-lfi\{B-^QyB-\ 

the m-th order Taylor expansion of l is 

m 

Tra{B + 0)-i = + Y,{-iy{B-'QyB-\ (16) 

It means that if we define Ui — Q [i — \,. . . ,to) then our approximation (Equation (|15p ) gives back 
the Taylor expansion. However, these ai parameters guarantee that our m-th order approximation 
gives the exact result if m > rankQ, while the Taylor expansion gives just an approximation for every 

TO. 

Remark 2.1. Assume that g is a metric on V, that is a bilinear symmetric map g :V xV which 
is non-degenarate: for every ^ v €V there exists a vector u such that g{u,v) ^ 0. Then for 
every v € V we have g{v, ■) G V* , and we can define an isomorphism 

g:V^V* v^g{v,-). 

Assume that A : V ^ V* is a. linear map and define A = g~^A which is a,V ^ V linear map. If A is 
invertible and perturbated to A' = A + W, then the above-mentioned Theorem and Corollary gives 
approximations (^')m^ fo^' (^0^^ ^^'^ '^'^ have the approximations {A')^g^^ for {A')^^. 
Remark 2.2. If dim!/ e {2, . . . , 10} and the matrix _B is a random matrix, fc e {2, . . . , 15} and the 
vectors (wi)i=i,...,fc and covectors (pi)i=i,...,fc are random vectors, and V is endowed with the Euclidean 
metric, then numerical simulations show, that the convergence of the given approximation is faster 
than the convergence of the Taylor expansion. We conjecture that this numerical observation is true 
in general settings too. 
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